Alfven waves as a driving mechanism in stellar winds 
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Alfven waves have been invoked as an important mechanism of particle acceleration in 
stellar winds of cool stars. After their identification in the solar wind they started to be 
studied in winds of stars located in different regions of the HR diagram. We discuss here 
Oh 1 some characteristics of these waves and we present a direct application in the acceleration 
q of late-type stellar winds. 

1. Introduction 

In 1942 combining fluid hydrodynamics and electromagnetism Hannes Alfven discovered 
a new mode of waves, thus opening a new field of research. He studied the mutual 
interaction between conducting fluid motion and electromagnetic fields and discovered 
these waves (Alfven, 1942) that later on were named Alfven waves. However, to be 
verified the waves needed to be detected in laboratory experiments. This was done in the 
late 40s and beginning of the 50s, in experiments with mercury and liquid sodium. It was 
discovered that in liquid metals the waves are strongly damped. In the 60s, experiments in 
plasma demonstrated weakly damped Alfven waves with the properties predicted by the 
theory. Since then this field has expanded mainly due to space studies and thermonuclear 
research, making possible, among other things, to create high temperature plasma in 
laboratories. 

An Alfven wave propagating in a plasma is a traveling oscillation of the ions and the 
magnetic field. The ion mass density provides the inertia and the magnetic field line 
tension provides the restoring force. The wave vector can either propagate in the parallel 
direction of the magnetic field or at oblique incidence. The waves efficiently carry energy 
and momentum along the magnetic field lines. 

Identification of Alfven waves in the solar wind by means of spacecraft measurements 
was first achieved in late 60s. In 2007, Tomczyk et al. (2007) reported the detection of 
Alfven waves in images of the solar corona with the Coronal Multi-Channel Polarimeter 
instrument at the National Solar Observatory, New Mexico. Also in 2007, several research 
groups reconfirmed that Alfven waves have sufficient energy to heat the solar corona and 
momentum to accelerate the solar wind (De Pontieu et al., 2007; Okamoto et al., 2007). 
The Sun is a benchmark for stellar research providing an example of a cool main sequence 
star losing mass. 

Since the discovery of the high velocity winds from O and B supergiant stars by Morton 
(1967) and his co-workers (Morton, Jenkins & Brooks, 1969), the theories of stellar winds 
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had their interest amplified. On another part of the HR diagram, the winds from evolved 
K and M stars, characterized by low velocities, demanded a different theoretical approach. 
Although presenting very different terminal velocities, in the range of 600 - 3500 km s _1 
for early-type and 10 - 100 km s" 1 for late-type K and M stars, these stars present a 
similarity: both luminous early and late-type stars are losing mass at a rate as high as 
10~ 5 M yr" 1 (Lamers & Cassinelli, 1999). 

The idea that a flux of Alfven waves could accelerate winds of stars in many regions 
of the HR diagram remount from a few decades ago (Hollweg, 1973, 1987; Hartmann 
& MacGregor, 1980; Holzer, Fla & Leer, 1983; Underhill 1983; Hartmann & Avrett, 
1984; Jatenco-Pereira and Opher, 1989a,b,c; dos Santos, Jatenco-Pereira and Opher, 1993; 
Jatenco-Pereira, Opher and Yamamoto, 1994; among others). The belief that Alfven 
waves are present in the magnetized winds of cool giant stars is supported by the fact 
that these waves are observed in the solar wind (Smith et al., 1995; Balogh et al., 1995). 
If there are oscillations in the magnetic field at the base of the wind, Alfven waves will be 
generated. As they propagate outward, the dissipation of their energy and the transfer 
of their momentum to the plasma can heat and accelerate the wind. As an example of 
how these waves can be used as an acceleration mechanism of a stellar wind, we show 
later on this paper a simple model to explain the mass-loss in a typical supergiant K5 
star (Section 3). In Section 2, we present the issue of the damping of the Alfven waves. 
In Section 4, we apply this model to a typical K5 supergiant star and in Section 5, we 
draw our conclusions. 

2. The damping mechanism for Alfven waves 

For cool red giant stars the flux of Alfven waves can provide an efficient driving mech- 
anism for the wind. In this mechanism it is the wave magnetic pressure that accelerates 
the wind. Hartmann & MacGregor (1980) showed that if the waves are not damped, the 
resultant terminal velocities are much higher than observed. Assuming that the waves are 
damped with a constant damping length, they obtained terminal velocities and mass-loss 
rate consistent with the observations. 

In order to improve the model presented by Hartmann & MacGregor (1980), Jatenco- 
Pereira & Opher (1989a) discarded the use of an artificial constant damping length for 
the waves and considered three different physical damping mechanisms: (i) the nonlinear 
damping, (ii) the resonant surface damping, and (iii) the turbulent damping. They showed 
that, for an isothermal atmosphere with T ~ 10 4 K, their model can reproduce the 
observed large mass-loss rates and the small ratios between the terminal velocity and the 
escape velocity at the base (■UooA'eo) °f these stellar winds. 

To illustrate this mechanism, we present a simplified version of this model (see Vidotto 
& Jatenco-Pereira, 2006; Falceta-Gongalves, Vidotto & Jatenco-Pereira, 2006 for a more 
complete analysis). We focus our study on the resonant surface damping for the Alfven 
waves. A surface wave exists at the interface separating two environments, e.g., open 
magnetic flux tubes. Magnetohydrodynamic surface waves may decay through a process 
called "resonance absorption" (Hollweg, 1987), that leads to a concentration of the surface 
wave energy into the thin resonant layer. In this layer, the energy may be dissipated locally 
as heat or concentrated in the form of turbulent Alfven waves which may propagate further 
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and be dissipated elsewhere, or concentrated in other forms, such as high energy particles, 
electric currents, etc. In the present paper we suggest that resonant absorption dissipates 
the energy locally as heat with the damping length L given by (Jatenco-Pereira & Opher, 
1989a) 



where Lq is the damping length of the wave at the base of the wind, M — u/va is the 
Alfven-Mach number, the Alfven speed is given by va = (B/y/Airp) and A(r) is the cross- 
section of a flux tube at a radial distance r. The index "0" indicates that the parameter 
is being evaluated at the stellar surface, i.e., at r = r . If a large damping length is 
considered (i.e., there is a smooth damping of the wave throughout the wind), there will 
be a non negligible outward force acting beyond the critical point of the wind, creating 
thus a high-speed wind. Hence, the damping length cannot be too large. This will also 
assure that the WKB approximation is valid in the wind we study. 



The first attempts in accelerating winds by means of Alfven waves considered that 
the waves were propagating in an isothermal ambient. This assumption simplifies the 
equations, since the energy equation simply reads T = constant. However, if we are 
interested in modeling the temperature profile of the wind, there is a need to solve a 
consistent energy equation simultaneously with the mass and momentum equations of 
the wind. In this way, we are able to obtain both the temperature and the velocity 
profiles of the wind and compare them with observations. 

3.1. Geometry of the Magnetic Field Lines 

For a long time, it is known that the coronal holes are the source of the high-speed 
solar wind streams at the Earth's orbit (e.g., Peter & Judge, 1999). Solar coronal holes 
observations show a superradial expansion of the magnetic flux tubes: the area expansion 
of a solar coronal hole is 5 — 13 times greater than for a radial expansion (Lie-Svendsen, 
Hansteen & Leer, 2002; Esser et al., 2005; Tu et al., 2005). Jatenco-Pereira & Opher 
(1989a) implemented the Alfven wave driven wind model by taking into account our 
knowledge on coronal holes in the Sun. We apply the same assumption for cool giant 
stars. For simplicity, we assume that in the wind of a typical K5 supergiant star, the 
non-radial expansion factor is 



where Q and f2 are the solid angles measured at the stellar surface and at the edge of 
the coronal hole, respectively. Hence, in order to model the diverging geometry for the 
magnetic field, we follow the idea given by Kuin and Hearn (1982). Thus, the cross-section 
of a flux tube at a radial distance r is given by 




(1) 



3. Model 
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Figure 1. Magnetic field geometry used in the model (not in scale). The stellar radius is 
denoted by r$. The magnetic field geometry is divergent up to the transition radius r t , 
i.e., it has a superradial index S > 2. Beyond r t , the magnetic field lines become radial 
and the index equals 2. 



where S, the superradial index, is a parameter that determines the divergence of the 
geometry from the stellar radius r up to the transition radius r t , as shown in Figure [TJ 
For a given S > 2, r t is obtained by 



Q A(r t )/r 2 fr t \ s ~ 2 



Q A{r Q )/rl Vr 
Hence, we have 

rt = F l/{S-2) r ^ 



(3) 



(4) 



Considering this geometry, the conservation of magnetic flux yields a magnetic field in- 
tensity of 



B(r) = B< 



A(r ) 
' A(r) ' 



(5) 



A detailed study of the diverging geometry of the magnetic field in the winds of late-type 
stars can be found in Jatenco-Pereira & Opher (1989a). 

3.2. Wind Equations 

The equation of mass continuity expresses the conservation of mass. If the flow velocity 
is denoted by u and the gas density is p, we have in steady-state 

puA(r)=C. (6) 
The constant C in equation §6§ is computed at the wind base, i.e., C = p u A(r ). 
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Assuming a steady flow, the equation of motion is written as 
du GM+ dP d ({SB) 2 ) 



pu— = -p — — 

dr r z 



dr dr 8n 



(7) 



where —pGM±/r 2 is the gravitational attraction force, —dP/dr is the pressure gradient, 
and d[((5 B) 2 ) / (87r)]dr is the Alfven wave magnetic pressure gradient. The parameter SB 
is the magnetic field amplitude of the wave. It is related to the energy density e of the 
wave by 



e = 



((SB)'' 



4n 



(8) 



The energy equation is determined through the balance between wave heating, adiabatic 
expansion, and radiative cooling (Hartmann, Edwards, & Avrett, 1982). Assuming an 
ideal gas, we can write the gas pressure as P = phsT/m, where ks is the Boltzmann 
constant, T is the gas temperature, and m is the mean mass per particle. Thus, neglecting 
conduction, we write the energy equation as 



pu 



d_ 

dr 



5 k B T GM* 



u 

T + 2 in 



ude , „ „ . 



0) 



The term (u/2)de/dr is the rate at which the waves do work on the gas. Q is the wave 
heating rate, i.e., the rate at which the gas is being heated due to dissipation of wave 
energy, and Pr is the radiative cooling rate, both in erg cm~ 3 s -1 . The wave heating can 
be written as (Hollweg, 1973) 



Q 



(u + V A ) 



and the radiative cooling is given by 
P R = A n e n H , 



(10) 



11^ 



where n e is the electron density, % is the hydrogen density and A is the radiative loss 
function. Here, we adopt A given by Schmutzler & Tscharnuter (1993) and calculate n e 
as Hartmann & MacGregor (1980). 

If the waves are damped, the wave energy is dissipated and we write 



e = e - 



Mn /1 + Mn 



M Vl + M 



exp 



(12) 



The initial energy density e and the initial wave flux are related to each other by 
(Jatenco-Pereira & Opher, 1989a) 



<pA = toVAo (l + \ M o) ■ 

From the equations above, we write the temperature variation as 
r(Q - Pr) 



dT 

dr 



IT 

3 r 



pu(ksT /m) 



' z+ r_du\ 
u dr I 



(13) 



(14) 
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where, we defined Z as 
Z = 



S if r < Tt 
2 if r > rt . 



As e = p((5v) 2 ), where Sv is the velocity fluctuation of the wave, the velocity gradient is 
given by 



1 du 
u dr 



9 5 k B T 
u 



3 m 



n + 3M\~ 


z 


V 1 + M )_ 


r 



5k B T 2r(Q-P R ) GM, 



3 m 3 Zpu rZ 

((5v) 2 } ((H 2 )/1 + 3M 
-r + 



2LZ 4 V 1 + M 



(15) 



Thus, in order to obtain self-consistently the velocity and the temperature profiles of 
the wind, it is necessary to solve equations (jSJ), (THj) . and (|T5|) simultaneously. 

4. Results and discussions 

The model is applied here to the wind of a typical K5 supergiant star with a mass of 
16 M and a radius of 400 R . We adopt a temperature of 3500 K at the base of the wind, 
and assume that the magnetic field intensity is ~ 10 G at this location. In our model, 
we adopt an initial wave flux of <pA ~ 10 6 erg cm -2 s _1 , which is of the same order of 
magnitude as the one estimated for the Sun (Banerjee et al., 1998), and take L = 0.2 r 
and S = 5.2. The initial density is po — 1-07 x 10~ 13 g cm -3 . 

We perform the calculations until r = 300 r , where the terminal velocity has been 
already reached. 

4.1. Typical K5 supergiant star 

Similarly to the wind first studied by Parker (1960), the momentum equation (eq. [T5] ) 
has a critical point when the term in brackets on the left-hand side becomes zero, and hence 
the bracketed term on the right-hand side must necessarily be zero too. The requirement 
that the wind velocity should increase through the critical point determines the initial 
velocity at the base of the wind, i.e., du(r = ro)/dr > 0. Hence, to find the critical solution 
for a given set of initial conditions, we iterate the initial velocity until the solution passes 
through the critical point. Together with the initial density, the initial velocity determines 
the mass- loss rate. 

Our results are shown in Figure 2, where the velocity and the temperature profiles 
obtained were plotted for r < 30 tq. We obtained a terminal velocity of Uoo ~ 63 km s _1 ~ 
0.5f eo and a mass-loss rate of M ~ 1.2 x 10~ 7 M Q yr _1 . These results compares favorably 
with observations: observationally, the terminal velocity should be lower than the surface 
escape velocity in a ratio around 1:2, and the mass-loss rate should range from 10 -8 to 
lO" 5 M yr" 1 (Linsky 1998). 

According to Roddier & Roddier (1985), the CII emission detected in cool winds of giant 
stars is associated to regions with temperatures of 7000 — 9000 K. Carpenter, Robinson, 
& Judge (1995) obtained, for the wind of 7 Cru (M3.4III), a maximum temperature of 
~ 9000 K. Carpenter et al. (1999) found a maximum temperature of ~ 10000 K for the 
wind of A Vel (K5Ib). For ( Aur (K4Ib), Eaton (1993) found that the wind temperature 
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Figure 2. Velocity (solid line) and temperature (dashed line) profiles obtained for the 
wind of a typical K5 supergiant star. 



reaches 12000 K at ~ 2r . We can see in Figure 2 an abrupt rise in the temperature of 
the wind: the temperature rises from 3.5 x 10 3 K to ~ 10 4 K in the range from 1 r to 
~ 2 r . This heating is mainly due to the dissipation of the Alfven waves, ie, caused by 
the wave damping. At large distances, when almost all the flux of Alfven waves has been 
dissipated, the flow expands adiabatically. Hence, in the absence of strong heating, the 
temperature tends to fall monotonically with the adiabatic exponent 4/3. 

As to exemplify the effects of the choice of L in our model, we have run a second model 
taking into account a different value of initial damping length: L = 0.4 r . Comparing 
this second model with the previous one, we obtained a value of ~ 100 km/s while for 
the first model a flow with a terminal velocity of ~ 63 km/s was achieved. This situation 
is explained by the position where a large fraction of the energy is deposited in the wind. 
If energy is mainly deposited in the super-critical part of the flow, as is the case for 
the second model, the initial structure of the wind is not modified, i.e., it maintains the 
same mass-loss rate, but the super-critical part is accelerated. Comparing the mass-loss 
rate, we obtained for the first model M ~ 1.2 x 10~ 7 M Q yr _1 , while for the second 
M ~ 1.3 X 10" 7 M Q yr" 1 . 

5. Conclusions 

We presented here some characteristics of Alfven waves and showed how these waves 
can be used as an acceleration mechanism of stellar winds. As a particular case, we 
evaluated the wind temperature and velocity profiles of a typical K5 supergiant star using 
an outward-directed flux of Alfven waves as the main acceleration mechanism of the wind. 
For a typical K5 supergiant star, we obtained a high mass-loss rate and a low terminal 
velocity consistent with observations. 
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It is important to notice that after the recent report on the unambiguously detection of 
Alfven waves in the solar corona, these waves will increase their applicability in the Sun 
and in several astrophysical environments. 
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